Paired guide RNA CRISPR-Cas9 screening for protein-coding genes and lncRNAs involved in transdifferentiation of human B-cells to macrophages.

CRISPR-Cas9 screening libraries have arisen as a powerful tool to identify protein-coding (pc) and non-coding genes playing a role along different processes. In particular, the usage of a nuclease active Cas9 coupled to a single gRNA has proven to efficiently impair the expression of pc-genes by generating deleterious frameshifts. Here, we first demonstrate that targeting the same gene simultaneously with two guide RNAs (paired guide RNAs, pgRNAs) synergistically enhances the capacity of the CRISPR-Cas9 system to knock out pc-genes. We next design a library to target, in parallel, pc-genes and lncRNAs known to change expression during the transdifferentiation from pre-B cells to macrophages. We show that this system is able to identify known players in this process, and also predicts 26 potential novel ones, of which we select four (two pc-genes and two lncRNAs) for deeper characterization. Our results suggest that in the case of the candidate lncRNAs, their impact in transdifferentiation may be actually mediated by enhancer regions at the targeted loci, rather than by the lncRNA transcripts themselves. The CRISPR-Cas9 coupled to a pgRNAs system is, therefore, a suitable tool to simultaneously target pc-genes and lncRNAs for genomic perturbation assays.

Roadmap for the use of base editors to decipher drug mechanism of action.

CRISPR base editors are powerful tools for large-scale mutagenesis studies. This kind of approach can elucidate the mechanism of action of compounds, a key process in drug discovery. Here, we explore the utility of base editors in an early drug discovery context focusing on G-protein coupled receptors. A pooled mutagenesis screening framework was set up based on a modified version of the CRISPR-X base editor system. We determine optimized experimental conditions for mutagenesis where sgRNAs are delivered by cell transfection or viral infection over extended time periods (>14 days), resulting in high mutagenesis produced in a short region located at -4/+8 nucleotides with respect to the sgRNA match. The ß2 Adrenergic Receptor (B2AR) was targeted in this way employing a 6xCRE-mCherry reporter system to monitor its response to isoproterenol. The results of our screening indicate that residue 184 of B2AR is crucial for its activation. Based on our experience, we outline the crucial points to consider when designing and performing CRISPR-based pooled mutagenesis screening, including the typical technical hurdles encountered when studying compound pharmacology.

3'Pool-seq: an optimized cost-efficient and scalable method of whole-transcriptome gene expression profiling.

BACKGROUND: The advent of Next Generation Sequencing has allowed transcriptomes to be profiled with unprecedented accuracy, but the high costs of full-length mRNA sequencing have posed a limit on the accessibility and scalability of the technology. To address this, we developed 3'Pool-seq: a simple, cost-effective, and scalable RNA-seq method that focuses sequencing to the 3'-end of mRNA. We drew from aspects of SMART-seq, Drop-seq, and TruSeq to implement an easy workflow, and optimized parameters such as input RNA concentrations, tagmentation conditions, and read depth specifically for bulk-RNA. RESULTS: Thorough optimization resulted in a protocol that takes less than 12 h to perform, does not require custom sequencing primers or instrumentation, and cuts over 90% of the costs associated with TruSeq, while still achieving accurate gene expression quantification (Pearson's correlation coefficient with ERCC theoretical concentration r = 0.96) and differential gene detection (ROC analysis of 3'Pool-seq compared to TruSeq AUC = 0.921). The 3'Pool-seq dual indexing scheme was further adapted for a 96-well plate format, and ERCC spike-ins were used to correct for potential row or column pooling effects. Transcriptional profiling of troglitazone and pioglitazone treatments at multiple doses and time points in HepG2 cells was then used to show how 3'Pool-seq could distinguish the two molecules based on their molecular signatures. CONCLUSIONS: 3'Pool-seq can accurately detect gene expression at a level that is on par with TruSeq, at one tenth of the total cost. Furthermore, its unprecedented TruSeq/Nextera hybrid indexing scheme and streamlined workflow can be applied in several different formats, including 96-well plates, which allows users to thoroughly evaluate biological systems under several conditions and timepoints. Care must be taken regarding experimental design and plate layout such that potential pooling effects can be accounted for and corrected. Lastly, further studies using multiple sets of ERCC spike-ins may be used to simulate differential gene expression in a system with known ground-state values.

SRSF1 and PTBP1 Are trans-Acting Factors That Suppress the Formation of a CD33 Splicing Isoform Linked to Alzheimer's Disease Risk.

A single nucleotide polymorphism (SNP) in exon 2 of the CD33 gene is associated with reduced susceptibility to late-onset Alzheimer's disease (AD) and causal for elevated mRNA lacking exon 2. In contrast to full-length CD33, transcripts lacking exon 2 result in CD33 protein unable to suppress activation responses in myeloid cells, including microglia. Currently, little is known about the regulation of CD33 exon 2 splicing. Using functional genomics and proteomic approaches, we found that SRSF1 and PTBP1 act as splicing enhancers to increase CD33 exon 2 inclusion in mRNA. Binding of PTBP1 to RNA sequences proximal to the intron 1-exon 2 splice junction is altered by the SNP and represents a potential mechanism behind the SNP-genotype dependent alternative splicing. Our studies also reveal that binding of SRSF1 to the CD33 RNA is not altered by the SNP genotype. Instead, a putative SRSF1 binding sequence at the 3' end of exon 2 directs CD33 exon 2 inclusion into the mRNA, indicating that PTBP1 and SRSF1 promote full-length isoform expression through different mechanisms. Our findings shed light on molecular interactions that regulate CD33 exon 2 splicing, ultimately impacting receptor expression on the cell surface. These data aid in the understanding of CD33's regulation of microglial signaling underpinning the AD genetic associations.

Scalable Design of Paired CRISPR Guide RNAs for Genomic Deletion.

CRISPR-Cas9 technology can be used to engineer precise genomic deletions with pairs of single guide RNAs (sgRNAs). This approach has been widely adopted for diverse applications, from disease modelling of individual loci, to parallelized loss-of-function screens of thousands of regulatory elements. However, no solution has been presented for the unique bioinformatic design requirements of CRISPR deletion. We here present CRISPETa, a pipeline for flexible and scalable paired sgRNA design based on an empirical scoring model. Multiple sgRNA pairs are returned for each target, and any number of targets can be analyzed in parallel, making CRISPETa equally useful for focussed or high-throughput studies. Fast run-times are achieved using a pre-computed off-target database. sgRNA pair designs are output in a convenient format for visualisation and oligonucleotide ordering. We present pre-designed, high-coverage library designs for entire classes of protein-coding and non-coding elements in human, mouse, zebrafish, Drosophila melanogaster and Caenorhabditis elegans. In human cells, we reproducibly observe deletion efficiencies of ≥50% for CRISPETa designs targeting an enhancer and exonic fragment of the MALAT1 oncogene. In the latter case, deletion results in production of desired, truncated RNA. CRISPETa will be useful for researchers seeking to harness CRISPR for targeted genomic deletion, in a variety of model organisms, from single-target to high-throughput scales.

Human selenoprotein P and S variant mRNAs with different numbers of SECIS elements and inferences from mutant mice of the roles of multiple SECIS elements.

Dynamic redefinition of the 10 UGAs in human and mouse selenoprotein P (Sepp1) mRNAs to specify selenocysteine instead of termination involves two 3' UTR structural elements (SECIS) and is regulated by selenium availability. In addition to the previously known human Sepp1 mRNA poly(A) addition site just 3' of SECIS 2, two further sites were identified with one resulting in 10-25% of the mRNA lacking SECIS 2. To address function, mutant mice were generated with either SECIS 1 or SECIS 2 deleted or with the first UGA substituted with a serine codon. They were fed on either high or selenium-deficient diets. The mutants had very different effects on the proportions of shorter and longer product Sepp1 protein isoforms isolated from plasma, and on viability. Spatially and functionally distinctive effects of the two SECIS elements on UGA decoding were inferred. We also bioinformatically identify two selenoprotein S mRNAs with different 5' sequences predicted to yield products with different N-termini. These results provide insights into SECIS function and mRNA processing in selenoprotein isoform diversity.

DECKO: Single-oligo, dual-CRISPR deletion of genomic elements including long non-coding RNAs.

BACKGROUND: CRISPR genome-editing technology makes it possible to quickly and cheaply delete non-protein-coding regulatory elements. We present a vector system adapted for this purpose called DECKO (Double Excision CRISPR Knockout), which applies a simple two-step cloning to generate lentiviral vectors expressing two guide RNAs (gRNAs) simultaneously. The key feature of DECKO is its use of a single 165 bp starting oligonucleotide carrying the variable sequences of both gRNAs, making it fully scalable from single-locus studies to complex library cloning. RESULTS: We apply DECKO to deleting the promoters of one protein-coding gene and two oncogenic lncRNAs, UCA1 and the highly-expressed MALAT1, focus of many previous studies employing RNA interference approaches. DECKO successfully deleted genomic fragments ranging in size from 100 to 3000 bp in four human cell lines. Using a clone-derivation workflow lasting approximately 20 days, we obtained 9 homozygous and 17 heterozygous promoter knockouts in three human cell lines. Frequent target region inversions were observed. These clones have reductions in steady-state MALAT1 RNA levels of up to 98 % and display reduced proliferation rates. CONCLUSIONS: We present a dual CRISPR tool, DECKO, which is cloned using a single starting oligonucleotide, thereby affording simplicity and scalability to CRISPR knockout studies of non-coding genomic elements, including long non-coding RNAs.

Splicing of a non-coding antisense transcript controls LEF1 gene expression.

In this report we have analyzed the role of antisense transcription in the control of LEF1 transcription factor expression. A natural antisense transcript (NAT) is transcribed from a promoter present in the first intron of LEF1 gene and undergoes splicing in mesenchymal cells. Although this locus is silent in epithelial cells, and neither NAT transcript nor LEF1 mRNA are expressed, in cell lines with an intermediate epithelial-mesenchymal phenotype presenting low LEF1 expression, the NAT is synthesized and remains unprocessed. Contrarily to the spliced NAT, this unspliced NAT down-regulates the main LEF1 promoter activity and attenuates LEF1 mRNA transcription. Unspliced LEF1 NAT interacts with LEF1 promoter and facilitates PRC2 binding to the LEF1 promoter and trimethylation of lysine 27 in histone 3. Expression of the spliced form of LEF1 NAT in trans prevents the action of unspliced NAT by competing for interaction with the promoter. Thus, these results indicate that LEF1 gene expression is attenuated by an antisense non-coding RNA and that this NAT function is regulated by the balance between its spliced and unspliced forms.

African signatures of recent positive selection in human FOXI1.

BACKGROUND: The human FOXI1 gene codes for a transcription factor involved in the physiology of the inner ear, testis, and kidney. Using three interspecies comparisons, it has been suggested that this may be a gene under human-specific selection. We sought to confirm this finding by using an extended set of orthologous sequences. Additionally, we explored for signals of natural selection within humans by sequencing the gene in 20 Europeans, 20 East Asians and 20 Yorubas and by analysing SNP variation in a 2 Mb region centered on FOXI1 in 39 worldwide human populations from the HGDP-CEPH diversity panel. RESULTS: The genome sequences recently available from other primate and non-primate species showed that FOXI1 divergence patterns are compatible with neutral evolution. Sequence-based neutrality tests were not significant in Europeans, East Asians or Yorubas. However, the Long Range Haplotype (LRH) test, as well as the iHS and XP-Rsb statistics revealed significantly extended tracks of homozygosity around FOXI1 in Africa, suggesting a recent episode of positive selection acting on this gene. A functionally relevant SNP, as well as several SNPs either on the putatively selected core haplotypes or with significant iHS or XP-Rsb values, displayed allele frequencies strongly correlated with the absolute geographical latitude of the populations sampled. CONCLUSIONS: We present evidence for recent positive selection in the FOXI1 gene region in Africa. Climate might be related to this recent adaptive event in humans. Of the multiple functions of FOXI1, its role in kidney-mediated water-electrolyte homeostasis is the most obvious candidate for explaining a climate-related adaptation.